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Two  stellar  components  in  the  halo 
of  the  Milky  Way 

Daniela  Carollo1’2’3,5,  Timothy  C.  Beers2,3,  Young  Sun  Lee2'3,  Masashi  Chiba4,  John  E.  Norris5,  Ronald  Wilhelm6, 
Thirupathi  Sivarani2’3,  Brian  Marsteller2,3,  Jeffrey  A.  Munn7,  Coryn  A.  L.  Bailer-Jones8,  Paola  Re  Fiorentin8,9 
&  Donald  G.  York10,11 

The  halo  of  the  Milky  Way  provides  unique  elemental  abundance  and  kinematic  information  on  the  first  objects  to  form  in 
the  Universe,  and  this  information  can  be  used  to  tightly  constrain  models  of  galaxy  formation  and  evolution.  Although 
the  halo  was  once  considered  a  single  component,  evidence  for  its  dichotomy  has  slowly  emerged  in  recent  years  from 
inspection  of  small  samples  of  halo  objects.  Here  we  show  that  the  halo  is  indeed  clearly  divisible  into  two  broadly 
overlapping  structural  components— an  inner  and  an  outer  halo— that  exhibit  different  spatial  density  profiles,  stellar  orbits 
and  stellar  metallicities  (abundances  of  elements  heavier  than  helium).  The  inner  halo  has  a  modest  net  prograde  rotation, 
whereas  the  outer  halo  exhibits  a  net  retrograde  rotation  and  a  peak  metallicity  one-third  that  of  the  inner  halo.  These 
properties  indicate  that  the  individual  halo  components  probably  formed  in  fundamentally  different  ways,  through 
successive  dissipational  (inner)  and  dissipationless  (outer)  mergers  and  tidal  disruption  of  proto-Galactic  clumps. 


Astronomers  have  long  sought  to  constrain  models  for  the  formation 
and  evolution  of  the  Milky  Way  (our  Galaxy)  on  the  basis  of  observa¬ 
tions  of  the  stellar  and  globular  cluster  populations  that  it  contains. 
These  populations  are  traditionally  defined  as  samples  of  objects  that 
exhibit  common  spatial  distributions,  kinematics  and  metallicities 
(the  age  of  a  population,  when  available,  is  also  sometimes  used). 
Metallicity  is  taken  by  astronomers  to  represent  the  abundances  of 
elements  heavier  than  helium,  which  are  only  created  by  nucleo¬ 
synthesis  in  stars — either  internally  via  nuclear  burning  in  their 
cores  or  externally  during  explosive  nucleosynthesis  at  the  end  of 
their  lives.  The  earliest  generations  of  stars  have  the  lowest  metal- 
licities,  because  the  gas  from  which  they  formed  had  not  been 
enriched  in  heavy  elements  created  by  previous  stars  and  distributed 
throughout  the  primordial  interstellar  medium  by  stellar  winds  and 
supernovae. 

Previous  work  has  provided  evidence  that  the  halo  of  the  Milky 
Way  may  not  comprise  a  single  population,  primarily  from  analysis 
of  the  spatial  profiles  (or  inferred  spatial  profiles)  of  halo  objects1-4.  A 
recent  example  of  such  an  analysis  is  the  observation  of  two  different 
spatial  density  profiles  for  distinct  classes  of  RR  Lyrae  variable  stars  in 
the  halo5.  In  addition,  tentative  claims  for  a  net  retrograde  motion  of 
halo  objects  by  previous  authors  supports  the  existence  of  a  likely 
dual-component  halo6-10.  The  central  difficulty  in  establishing  with 
confidence  whether  or  not  a  dichotomy  of  the  halo  populations  exists 
is  that  the  past  samples  of  tracer  objects  have  been  quite  small,  and 
usually  suitable  only  for  consideration  of  a  limited  number  of  the 
expected  signatures  of  its  presence. 

In  the  present  work,  we  examine  this  question  in  detail  using  a 
large,  homogeneously  selected  and  analysed  sample  of  over  20,000 
stars,  originally  obtained  as  calibration  data  during  the  course  of  the 
Sloan  Digital  Sky  Survey  (SDSS)11.  Although  there  are  many  possible 
alternative  (and  more  complex)  models  that  might  be  considered, 


multiple  lines  of  evidence  derived  from  these  data  clearly  confirm 
that  the  halo  can  be  resolved  into  (at  least)  two  primary  populations, 
the  inner  and  the  outer  halo,  with  very  different  observed  properties. 

We  find  that  the  inner-halo  component  of  the  Milky  Way  domi¬ 
nates  the  population  of  halo  stars  found  at  distances  up  to  10-15  kpc 
from  the  Galactic  Centre  (including  the  solar  neighbourhood).  An 
outer-halo  component  dominates  in  the  regions  beyond  15-20  kpc. 
We  show  the  inner  halo  to  be  a  population  of  stars  that  are  non- 
spherically  distributed  about  the  centre  of  the  Galaxy,  with  an  inferred 
axial  ratio  of  the  order  of  ~0.6.  Inner-halo  stars  possess  generally  high 
orbital  eccentricities,  and  exhibit  a  modest  prograde  rotation 
(between  0  and  50kms  ’)  around  the  centre  of  the  Galaxy  (see 
Supplementary  Table  1).  The  distribution  of  metallicities  for  stars  in 
the  inner  halo  peaks  at  [Fe/H]  =  — 1.6,  with  tails  extending  to  higher 
and  lower  metallicities.  (Here  metallicity  is  defined  as  [A/B]  = 
log10(NA/NB)  —  oglo(NA/NB)0>  where  NA  and  NB  represent  the 
number  density  of  atoms  of  elements  A  and  B,  and  the  subscript  © 
indicates  solar  values.)  The  outer  halo,  by  contrast,  comprises  stars 
that  exhibit  a  much  more  spherical  spatial  distribution,  with  an  axial 
ratio  of  —0.9-1 .0.  Outer-halo  stars  cover  a  wide  range  of  orbital 
eccentricities,  including  many  with  lower  eccentricity  orbits  than 
found  for  most  stars  associated  with  the  inner  halo,  and  exhibit  a  clear 
retrograde  net  rotation  (between  -40  and  -70kms  ’)  about  the 
centre  of  the  Galaxy.  The  metallicity  distribution  function  (MDF)  of 
the  outer  halo  peaks  at  lower  metallicity  than  that  of  the  inner 
halo,  around  [Fe/H]  =  —2.2,  and  includes  a  larger  fraction  of  low- 
metallicity  stars  than  does  the  MDF  of  the  inner-halo  population. 

Evidence  for  the  dichotomy  of  the  halo 

The  spectroscopy,  photometry  and  astrometry  for  our  large  sample 
of  stars  were  obtained  from  observations  carried  out  with  the  Apache 
Point  2.5-m  SDSS  telescope;  these  data  are  publicly  available  as  Data 
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Figure  1  |  The  spatial  distribution  of  the  stars  analysed  in  the  present 
sample.  The  distribution  of  the  full  sample  of  20,236  unique  SDSS  (Data 
Release  512)  spectrophotometric  and  telluric  calibration  stars  in  the  Z-R 
plane  is  shown,  where  Z  is  the  derived  distance  from  the  Galactic  plane  in  the 
vertical  direction  and  R  is  the  derived  distance  from  the  centre  of  the  Galaxy 
projected  onto  this  plane.  The  dashed  blue  line  represents  the  Galactic  plane, 
while  the  filled  orange  dot  is  the  position  of  the  Sun,  at  Z  =  0  kpc  and 
R  =  8.5  kpc.  The  ‘wedge  shape’  of  the  selection  area  is  the  result  of  limits  of 
the  SDSS  footprint  in  Galactic  latitude.  The  red  points  indicate  the  10,123 
stars  that  satisfy  our  criteria  for  a  local  sample  of  stars,  having 
7  kpc  <  R  <  10  kpc,  with  distance  estimates  from  the  Sun  d  <  4  kpc,  and  with 
viable  measurements  of  stellar  parameters  and  proper  motions. 

Release  512.  Details  concerning  the  selection  of  the  stars  and  the 
measurement  of  their  parameters  (temperature,  surface  gravity  and 
metallicity,  [Fe/H]),  as  well  as  the  methods  used  to  obtain  their 


estimated  distances,  proper  motions  and  derived  kinematics,  can 
be  found  in  the  Supplementary  Information. 

Our  20,236  programme  stars  explore  distances  up  to  20  kpc  from 
the  Sun,  but  we  can  only  obtain  useful  estimates  of  the  full  space 
motions  (as  described  in  the  Supplementary  Information)  for  the 
subset  of  10,123  stars  in  a  local  volume  (up  to  4  kpc  from  the  Sun; 
Fig.  1).  The  restriction  of  the  sample  to  the  region  of  the  solar  neigh¬ 
bourhood  is  also  made  so  that  the  assumptions  going  into  the  kine¬ 
matic  calculations  are  best  satisfied.  Figure  2  shows  the  distribution 
of  [Fe/H]  for  different  cuts  in  the  V  velocity,  which  is  the  orbital 
component  that  measures  the  motion  of  a  star  (with  respect  to  the 
Local  Standard  of  Rest)  in  the  rotation  direction  of  the  Galaxy.  The 
transition  in  the  distribution  of  [Fe/H]  that  is  expected  as  one  sweeps 
from  stars  with  thick-disk-like  motions  to  stars  with  halo-like 
motions  is  clear.  However,  with  the  large  sample  of  stars  in  our 
sample,  it  is  possible  to  investigate  the  change  in  the  distribution  of 
[Fe/H]  for  stars  that  are  increasingly  more  retrograde,  as  well  as  for 
those  that  are  both  highly  retrograde  and  have  orbits  taking  them  to 
high  (the  maximum  distance  above  the  Galactic  plane  reached 
by  a  star  during  the  course  of  its  orbit  about  the  Galactic  Centre — the 
Supplementary  Information  describes  the  methods  used  to  derive 
this  fundamental  parameter).  This  figure  shows  that  stars  with  the 
most  retrograde  orbits,  and  those  that  reach  large  distances  in  their 
orbits  above  the  Galactic  plane,  exhibit  distributions  of  [Fe/H]  that 
peak  at  metallicities  between  —2.0  and  —2.2,  which  we  associate  with 
the  outer-halo  population.  The  inner-halo  population  dominates  the 
samples  of  stars  with  peak  metallicity  [Fe/H]  ~  — 1.6. 

Astronomers  have  long  debated  whether  there  might  exist  a  change 
in  the  rotational  properties  of  the  halo  of  the  Milky  Way  as  a  function 
of  distance  from  the  Galactic  Centre,  based  on  much  smaller  samples 
of  globular  clusters2,10  and  stars6-9,13,14  than  we  consider  here.  The 


[Fe/H] 


Figure  2  |  The  distribution  of  [Fe/H]  for  various  cuts  in  the  V  velocity  (in 
km  s-1),  the  component  of  orbital  motion  measured  with  respect  to  the 
Local  Standard  of  Rest.  The  Local  Standard  of  Rest  is  a  frame  in  which  the 
mean  space  motions  of  the  stars  in  the  solar  neighbourhood  average  to  zero. 
A  blue  dashed  line  at  [Fe/H]  =  —2.0  is  added  for  reference  in  all  three 
columns.  In  the  left-hand  column,  the  full  data  set  is  considered.  The  stars 
with  modestly  negative  V  velocities  in  the  upper  three  panels  are  dominated 
by  stars  from  the  thick-disk  (and  metal-weak  thick-disk)  populations,  with  a 
peak  metallicity  around  [Fe/H]  ~  —0.7.  Atransition  to  dominance  by  inner- 
and  outer-halo  population  stars  becomes  evident  for  V<  —  90  km  s  ;  in  the 
bottom  panel  of  this  column,  the  distribution  of  [Fe/H]  appears  similar  to 
what  in  the  past  was  considered  ‘the  halo’,  but  we  argue  results  from  a 


superposition  of  contributions  from  both  inner-  and  outer-halo 
populations.  In  the  middle  column,  the  large  numbers  of  stars  with 
V  <  —200  km  s  1  are  broken  into  smaller  ranges  in  V  velocity.  As  V 
becomes  increasingly  retrograde  ( V<  —220  km  s  1 ),  the  metallicity 
distribution  shifts  to  include  ever  larger  numbers  of  stars  with 
[Fe/H]  <  —2.0,  and  relatively  fewer  stars  with  [Fe/H]  ~  — 1.6.  The  same  V 
velocity  cuts  are  applied  in  the  right-hand  column,  but  only  for  stars  with 
Zmax  >  5  kpc,  in  order  to  decrease  the  contribution  from  inner-halo  stars. 
Although  fewer  stars  are  included,  the  increasing  dominance  of  stars  with 
[Fe/H]  <  —2.0  is  even  more  apparent.  We  associate  the  stars  with  the  most 
extreme  retrograde  orbits  (and  those  that  reach  far  above  the  Galactic  plane 
in  their  orbits)  with  the  outer-halo  population. 
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stellar  samples  were  obtained  with  selection  criteria  (for  example,  on 
the  basis  of  high  proper  motions  for  halo  stars  in  the  solar  neigh¬ 
bourhood7'13,14,  or  from  in  situ,  apparent-magnitude  limited  sur¬ 
veys6'9,15)  that  we  suggest  favoured  membership  in  one  or  the  other 
of  the  now  clearly  revealed  halo  components.  A  detailed  summary  of 
the  kinematics  of  our  programme  stars  is  presented  in  the  Supple¬ 
mentary  Information,  where  we  also  establish  consistency  between 
properties  obtained  through  techniques  based  on  full  space  motions 
and  those  based  on  radial  velocities  alone,  which  argues  against  the 
existence  of  any  large  systematic  errors  in  the  proper  motions. 

Table  1  summarizes  the  past  and  present  determinations  of  (  V0 ), 
the  mean  rotational  velocity  with  respect  to  the  Galactic  Centre, 
where  claims  for  a  retrograde  halo  have  been  made.  Previous  samples 
that  have  addressed  this  question  were  based  on  either  much  smaller 
total  numbers  of  objects  (with  the  limitation  that  they  could  not  well 
sample  both  the  inner-  and  outer-halo  populations),  did  not  have 
proper  motions  available  (or  only  highly  uncertain  ones),  or  were 
otherwise  restricted  due  to  the  selection  criteria  employed  (that  is, 
they  were  kinematically  biased16,  or  had  limited  sky  coverage,  ren¬ 
dering  them  sensitive  to  the  effects  of  individual  star  streams17).  The 
local  sample  of  SDSS  calibration  stars  we  have  assembled  does  not 
suffer  from  any  of  these  limitations.  The  retrograde  signatures  for 
stars  we  associate  with  the  outer-halo  population  are  robust  and 
highly  statistically  significant  (except  for  the  smallest  subsample). 
However,  even  our  precise  present  determination  of  the  net  retro¬ 
grade  rotation  of  the  outer  halo,  based  on  our  local  sample,  is  prob¬ 
ably  influenced  by  some  degree  of  overlap  between  outer-halo  stars 
with  those  from  the  inner-halo  population. 

The  distribution  of  [Fe/H]  for  stars  on  increasingly  retrograde 
orbits  about  the  Galactic  Centre  for  subsamples  that  reach  different 
distances  from  the  Galactic  plane  {Z^, ax)  is  shown  in  Fig.  3.  The  MDFs 
of  the  stars  with  Zmax  close  to  the  Galactic  plane  are  very  different 
from  those  whose  orbits  reach  farther  from  the  plane.  The  distri¬ 
bution  of  metallicity  clearly  shifts  to  lower  abundances  as  more  severe 
cuts  on  V0  or  Zmax  are  applied,  as  supported  by  rigorous  statistical 
tests.  We  conclude  that  the  halo  of  the  Galaxy  comprises  stars  with 
intrinsically  different  distributions  of  [Fe/H];  the  observed  changes 
in  the  MDF  of  halo  stars  with  V0  and  Zmax  would  not  be  expected 
if  the  halo  is  considered  as  a  single  entity.  The  Supplementary 


Information  presents  additional  observed  differences  in  the  ener¬ 
getics,  the  distribution  of  orbital  eccentricities,  and  changes  in  the 
nature  of  stellar  orbits  for  our  programme  stars  that  are  also  incon¬ 
sistent  with  a  single  halo  population. 

In  order  to  provide  confirmation  of  the  shift  in  the  MDF  inferred 
from  our  analysis  of  a  local  sample  of  stars,  we  also  examine  an 
auxiliary  sample  of  stars  that  are  at  present  located  much  farther 
from  the  Galactic  Centre.  This  sample  comprises  1,235  blue  hori¬ 
zontal-branch  stars  selected  from  the  SDSS18.  The  stars  cover  a  wide 
range  of  distances,  from  5  kpc  to  over  80  kpc  from  the  centre  of  the 
Galaxy.  Statistical  tests  strongly  reject  the  hypothesis  that  the  stars  at 
large  distances  from  the  Galactic  Centre  could  be  drawn  from  the 
same  parent  population  as  those  at  distances  close  to  the  Galactic 
Centre  (Fig.  4). 

It  has  been  shown,  on  the  basis  of  Jeans’  theorem19'20,  that  the 
global  structure  of  the  stellar  halo  can  be  recovered  from  local  kine¬ 
matic  information,  as  long  as  one  has  a  sufficiently  large  number  of 
stars  observed  in  the  solar  neighbourhood  that  explore  the  full  phase- 
space  distribution  of  the  pertinent  stellar  populations.  We  note  that 
the  actual  halo  systems  of  the  Galaxy  are  unlikely  to  be  in  well-mixed 
equilibrium  states.  However,  the  relaxation  process  is  very  slow  com¬ 
pared  to  the  orbital  periods  of  typical  stars,  so  Jeans’  theorem  and  the 
approach  based  on  it  remain  at  least  approximately  valid.  The  result 
of  this  exercise  for  our  large  sample  of  SDSS  calibration  stars,  over 
narrow  cuts  in  metallicity,  is  shown  in  Fig.  5.  The  observed  changes  in 
the  inferred  spatial  density  profiles  suggest  that  a  flattened  inner-halo 
population  dominates  locally  for  stars  with  [Fe/H]  >  —  2,  whereas 
the  outer-halo  population  has  a  nearly  spherical  distribution,  and 
dominates  at  distances  beyond  r~  15-20  kpc  (where  r  represents 
the  distance  from  the  Galactic  Centre),  as  well  as  locally  for  stars  with 
[Fe/H]  <  —2.0.  Variations  in  the  halo  spatial  profile  with  distance 
have  been  recognized  by  a  number  of  previous  authors1^1’15'20,  based 
on  samples  of  stars  that  are  one  to  two  orders  of  magnitude  smaller 
than  our  present  data  set. 

Implications  of  the  dichotomy  of  the  halo 

An  early  model  for  the  formation  of  the  Milky  Way,  based  on  the 
rapid  (a  few  hundred  million  years)  monolithic  collapse  of  a  gaseous 
proto-Galaxy21,  has  yielded  to  the  more  recent  idea  that  the  halo  of 


Table  1  |  Studies  claiming  a  retrograde  outer  halo 


Sample  and  selection  criteria 

N 

Additional 

0A> 

Method 

Source 

restrictions 

(kms  1) 

Globular  clusters  (non-kinematic) 

19 

'Young  halo' 

-64  ±74 

F&W 

Ref.  2 

Globular  clusters  (non-kinematic) 

20 

'Young  halo' 

-42  ±80 

F&W 

Ref.  10 

RR  Lyrae  stars  (non-kinematic) 

26 

Z  <  8  kpc 

-95  ±29 

FSM 

Ref.  9 

Field  subdwarfs  (kinematic) 

30 

Zmax  '>  5  kpc 

-45  ±22 

FSM 

Ref.  7 

Bias  corrected 

±24  ±13 

Ref.  16 

Field  horizontal-branch  stars  (non-kinematic) 

90 

[Fe/H]  <-1.6 

-93  ±36 

F&W 

Ref.  8 

Z  >4kpc 

Field  subdwarfs  (kinematic) 

101 

V<  -100  km  s_1 

-32  ±  10 

FSM 

Ref.  13 

[Fe/H]  <-1.8 

Field  F,  G,  K  dwarfs  (non-kinematic) 

250 

Z  >  5  kpc 

-55  ±  16 

FSM 

Ref.  6 

Field  F,  G  turnoff  (non-kinematic) 

2,228 

Zmax  >5  kpc 

-11  ±2 

FSM 

This  work 

200 

[Fe/H]  <-1.0 
[Fe/H]  <-2.2 

-41  ±  11 

771 

Zmax  >  10  l<PC 

-38  ±5 

94 

[Fe/H]  <-1.0 
[Fe/H]  <-2.2 

-71  ±  17 

371 

Zmax  >15  kpc 

-56  ±8 

54 

[Fe/H]  <-1.0 
[Fe/H]  <-2.2 

-71  ±25 

Previous  and  current  determinations  of  the  mean  rotational  velocity,  (V#),  and  the  error  in  the  mean  (a/ VN,  where  a  is  the  standard  deviation  and  N  is  the  number  of  stars),  for  samples  in  which  a 
counter-rotating  halo  has  been  claimed,  ordered  by  sample  size.  The  samples  listed  in  the  first  column  are  classified  as  to  whether  they  were  chosen  on  the  basis  of  high  proper  motions  (kinematic)  or 
not  (non-kinematic).  Restrictions  placed  on  each  sample  by  the  authors  are  listed  in  the  third  column  (see  original  papers  for  details).  The  method  of  analysis  used  for  each  determination  is  listed: 
F&W,  estimate  based  on  the  technique  of  Frenk  and  White48,  which  considers  distances  and  radial  velocities  alone,  under  the  assumption  of  a  cylindrically  symmetric  Galaxy;  FSM,  estimate  based  on 
consideration  of  the  full  space  motions,  which  requires  the  use  of  proper  motions,  as  well  as  distances  and  radial  velocities.  The  samples  analysed  with  the  F&W  approach  are  either  not  statistically 
different  from  zero  (refs  2, 10),  or  are  only  marginally  so  (2.6<r;  ref.  8).  Previous  samples  based  on  analysis  of  the  full  space  motions  vary  from  statistically  insignificant  (ref.  7),  to  just  over  30- 
significance  (refs  6,  9, 13),  owing  to  the  small  numbers  of  stars  considered.  Note  that  after  application  of  (uncertain)  corrections  for  kinematic  bias  (ref.  16),  the  retrograde  result  reported  in  ref.  7 
disappears  entirely.  One  can  assume  that  a  similar  outcome  might  apply  to  the  ref.  13  determination.  The  samples  of  ref.  6  and  ref.  9  are  both  selected  over  a  restricted  region  of  sky  (towards  the 
North  Galactic  Pole)  and  are  therefore  subject  to  possible  contamination  by  individual  stellar  streams.  All  but  two  subsamples  of  the  SDSS  calibration  star  sample  have  retrograde  signals  that  are 
significant  at  more  than  the  4<r  level.  The  subsample  at  Zmax  >  5  kpc  and  [Fe/H]  <  —1  is  likely  to  include  significant  contamination  from  inner-halo  stars. 
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Figure  3  |  The  distribution  of  [Fe/H]  for  the  stars  in  our  sample  on  highly 
retrograde  orbits.  Stars  from  the  disk  populations,  which  possess  prograde 
orbits,  cannot  be  present  in  this  plot.  The  panels  show  various  cuts  in  V0 ,  the 
rotational  velocity  with  respect  to  the  Galactic  Centre  in  a  cylindrical 
coordinate  system,  and  for  different  ranges  of  Zmax  (in  kpc).  A  blue  dashed 
line  at  [Fe/H]  =  —2.0  is  added  for  reference.  The  left-hand  column  applies 
for  stars  with  V0  <  — 100  km  s-1;  the  clearly  skewed  distribution  of  [Fe/H] 
exhibits  an  increased  contribution  from  lower  metallicity  stars  as  one 
progresses  from  the  low  ( Zmax  <  5  kpc)  to  the  high  (Zmax  >15  kpc) 
subsamples.  Simultaneously,  the  predominance  of  stars  from  the  inner-halo 
population,  with  peak  metallicity  at  [Fe/H]  ~  — 1.6,  decreases  in  relative 
strength,  and  shifts  to  lower  [Fe/H].  Similar  behaviours  are  seen  in  the 


middle  and  right-hand  columns,  which  correspond  to  cuts  on 
V0  <  — 150  km  s-1  and  —200  km  s-1,  respectively.  A  Kolmogorov-Smirnoff 
test  of  the  null  hypothesis  that  the  MDFs  of  stars  shown  in  the  lower  panels 
for  the  individual  cuts  on  V0  could  be  drawn  from  the  MDFs  of  the  same 
parent  population  as  those  shown  in  the  upper  panels,  against  an  alternative 
that  the  stars  are  drawn  from  more  metal-poor  parent  MDFs,  is  rejected  at 
high  levels  of  statistical  significance.  For  V0  <  — 100  km  s-1,  one-sided 
probabilities  less  than  0.0001  are  obtained  for  the  cuts  on  Zmax  >  5,  10  and 
15kpc.For  V0<  — 150  km  s-1,  one-sided  probabilities  of 0.0004, 0.0001  and 
0.0003  are  obtained  for  Zmax  >  5,  10  and  15  kpc,  respectively.  For 
V0  <  —  200  km  s-1,  one-sided  probabilities  of  0.014,  0.010  and  0.033  are 
obtained,  for  Zmax  >  5,  10  and  15  kpc,  respectively. 


the  Galaxy  was  assembled,  over  the  span  of  several  billion  years,  from 
smaller  proto-Galactic  clumps22.  This  hierarchical  assembly  model 
has  received  close  attention  in  recent  years,  in  part  because  it  fits  well 
with  the  prevailing  theory  for  the  formation  and  evolution  of  struc¬ 
ture  in  the  Universe,  based  on  the  early  collapse  of  ‘mini-haloes’  of 
cold  dark  matter  (CDM)23,24.  Modern  numerical  simulations  for  the 


assembly  of  large  spirals  based  on  CDM  cosmogonies  predict  that  the 
stars  in  the  haloes  of  galaxies  like  the  Milky  Way  might  be  composed 
of  the  shredded  stellar  debris  of  numerous  dwarf-like  galaxies  that 
have  been  torn  apart  by  tidal  interactions  with  their  parent 
galaxy24-27.  Recent  quantitative  analysis  of  the  amount  of  structure 
visible  in  the  halo  of  the  Galaxy  from  SDSS28-31  imaging  provides 
compelling  additional  evidence32.  Others  have  argued  that  some 
combination  of  a  monolithic  collapse  and  a  hierarchical  assembly 
model  may  be  necessary  to  fully  explain  the  observed  data2,15,33. 

Within  the  context  of  the  CDM  model,  the  formation  of  the  inner 
halo  may  be  understood  in  the  following  manner.  Low-mass  sub- 
Galactic  fragments  are  formed  at  an  early  stage.  These  fragments 


Figure  4  |  A  sample  of  blue  horizontal-branch  stars  exploring  much  larger 
distances  from  the  Galactic  Centre  than  the  SDSS  calibration  stars.  The 

distribution  of  [Fe/H]  is  shown  for  various  cuts  on  the  distance  from  the 
Galactic  Centre,  r,  in  kpc.  The  nature  of  the  MDF  appears  to  shift  from  the 
upper  two  panels,  which  exhibit  the  character  of  a  mixture  of  inner-  and 
outer-halo  populations,  over  to  a  unimodal  distribution  in  the  third  panel, 
centred  on  [Fe/H]  ~  —2.0.  The  most  distant  blue  horizontal  branch  (BHB) 
stars  in  the  lowest  panel  also  exhibit  the  appearance  of  a  mixture  of  the  two 
populations,  possibly  due  to  the  inclusion  of  inner-halo  stars  on  highly 
eccentric  orbits  that  take  them  far  from  the  Galactic  Centre.  The  peak  around 
[Fe/H]  =  —3.0  seen  in  several  of  the  panels  is  an  artefact  arising  from  the 
limit  of  the  metallicity  grid  that  is  used  for  abundance  determinations  of  the 
BHB  stars.  A  Kolmogorov-Smirnoff  test  of  the  null  hypothesis  that  the 
MDFs  of  stars  shown  in  the  lower  panels  for  the  individual  cuts  on 
Galactocentric  distance  r  could  be  drawn  from  the  same  parent  population  as 
the  stars  shown  in  the  first  panel,  against  an  alternative  that  the  stars  are 
drawn  from  more  metal-poor  parent  MDFs,  is  rejected  at  high  levels  of 
statistical  significance  (one-sided  probabilities  of  0.0262, 0.0005  and  0.0243, 
respectively,  for  the  three  higher  cuts  on  r).  The  fraction  of  stars  with 
metallicities  [Fe/H]  <  —2.0  (primarily  outer-halo  stars)  grows  from  31%  for 
stars  with  5  <  r  <  10  kpc  to  46%  for  stars  at  larger  Galactocentric  distances. 
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rapidly  merge  into  several  (in  many  simulations,  two26'34)  more- 
massive  clumps,  which  themselves  eventually  dissipatively  merge 
(owing  to  the  presence  of  gas  that  has  yet  to  form  stars).  The  essen¬ 
tially  radial  merger  of  the  few  resulting  massive  clumps  gives  rise  to 
the  dominance  of  the  high-eccentricity  orbits  for  stars  that  we  assign 
here  to  membership  in  the  inner  halo.  Star  formation  within  these 
massive  clumps  (both  pre-  and  post-merger)  would  drive  the  mean 
metallicity  to  higher  abundances.  This  is  followed  by  a  stage  of  adia¬ 
batic  compression  (flattening)  of  the  inner  halo  component  owing  to 
the  growth  of  a  massive  disk,  along  with  the  continued  accretion  of 
gas  onto  the  Galaxy34,35. 

The  fact  that  the  outer-halo  component  of  the  Milky  Way  exhibits 
a  net  retrograde  rotation  (and  a  different  distribution  of  overall 
orbital  properties),  as  found  here,  clearly  indicates  that  the  formation 
of  the  outer  halo  is  distinct  from  that  of  both  the  inner-halo  and  disk 
components.  We  suggest,  as  others  have  before,  that  the  outer-halo 
component  formed,  not  through  a  dissipative,  angular-momentum- 
conserving  contraction,  but  rather  through  dissipationless  chaotic 
merging  of  smaller  subsystems  within  a  pre-existing  dark-matter 
halo.  These  subsystems  would  be  expected  to  be  of  much  lower  mass, 
and  subject  to  tidal  disruption  in  the  outer  part  of  a  dark-matter  halo, 
before  they  fall  farther  into  the  inner  part.  As  candidate  (surviving) 
counterparts  for  such  subsystems,  one  might  consider  the  low- 
luminosity  dwarf  spheroidal  galaxies  surrounding  the  Galaxy,  in  par¬ 
ticular  the  most  extreme  cases  recently  identified  from  the  SDSS36,37. 
Subsystems  of  lower  mass,  and  by  inference,  even  lower  metallicity, 
may  indeed  be  destroyed  so  effectively  that  none  (or  very  few)  have 
survived  to  the  present  day.  If  so,  the  outer-halo  population  may  be 
assembled  from  relatively  more  metal-poor  stars,  following  the 
luminosity-metallicity  relationship  for  Local  Group  dwarf  galaxies38. 


The  net  retrograde  rotation  of  the  outer  halo  may  be  understood  in 
the  context  of  the  higher  efficiency  of  phase  mixing  for  the  orbits  of 
stars  that  are  stripped  from  subsystems  on  prograde,  rather  than 
retrograde,  orbits39,40. 

The  clear  difference  in  the  MDFs  of  the  two  halo  populations  we 
identify  also  suggests  that  the  lowest  metallicity  stars  in  the  Galaxy 
may  be  associated  with  the  outer  halo,  which  can  be  exploited  for 
future  directed  searches.  It  is  noteworthy  that  the  hyper  metal-poor 
stars  HE  0107  —  524041  and  HE  1327  — 232642,  both  of  which  have 
[Fe/H]  <  —5.0,  as  well  as  the  recently  discovered  ultra  metal-poor 
star  HE  0557— 484043,  with  [Fe/H]  =  —4.8,  are  either  located  greater 
than  lOkpc  away  (HE  0107  —  5240,  HE  0557—4840)  or  have  space 
motions  that  carry  them  far  out  into  the  Galaxy  (HE  1327  —  2326; 
A.  Frebel,  personal  communication). 

In  addition,  efforts  to  determine  the  primordial  lithium  abun¬ 
dance  from  observations  of  the  most  metal-poor  stars44  may  have 
inadvertently  mixed  samples  from  the  inner-  and  outer-halo  popula¬ 
tions;  such  stars  could  have  formed  and  evolved  in  rather  different 
astrophysical  environments.  The  inner/outer  halo  dichotomy  may 
also  have  an  impact  on  the  expected  numbers  of  carbon-enhanced 
metal-poor  stars  as  a  function  of  declining  metallicity45,  and  as  a 
function  of  distance  from  the  Galactic  plane46,47. 

Much  remains  to  be  learned  as  the  database  of  low-metallicity  stars 
continues  to  expand,  in  particular  from  those  stars  that  are  found  in 
distant  in  situ  samples,  or  from  those  nearby  stars  with  available 
proper  motions  that  indicate  membership  of  the  outer-halo  popu¬ 
lation.  We  look  forward  to  the  next  dramatic  increase  in  the  numbers 
of  very  metal-poor  stars  that  will  come  from  the  ongoing  stellar 
samples  from  SDSS,  and  in  particular  from  SEGUE,  the  Sloan 
Extension  for  Galactic  Understanding  and  Exploration. 


-1 .2  <  [Fe/H]  <  -1 .0 
30 


S  15 


-1 .4  <  [Fe/H]  <  -1 .2 


Figure  5  |  Equidensity  contours  of  the  reconstructed  global  density 
distributions  for  stars  in  our  sample  with  various  metallicities.  The  global 
density  distributions  are  constructed  from  the  sum  of  the  probability  density 
of  an  orbit  at  each  location  in  the  Z—R  plane,  with  a  weighting  factor  being 
inversely  proportional  to  the  corresponding  density  at  the  currently 
observed  position  of  the  star19,20.  High-density  regions  are  indicated  by 
redder  colours,  while  low-density  regions  are  indicated  by  bluer  colours  (a 
linear  density  scale  is  used).  Within  each  metallicity  cut,  the  apparent 
flattening  of  the  inner  regions  slowly  goes  over  to  a  more  spherical  shape 


with  increasing  distance.  As  one  progresses  from  the  more  metal-rich 
([Fe/H]  —  —1.0)  to  the  most  metal-poor  ([Fe/H]  <  —2.2)  subsets  of  these 
data,  the  overall  nature  of  the  equidensity  contours  also  changes  from  highly 
flattened  (axial  ratios  of  ~0.6),  to  more  spherical  (axial  ratio  of  —0.9).  This 
suggests  that  the  inner-  and  outer-halo  components  are  broadly  overlapping 
in  space  and  in  metallicity — the  inner-halo  population  is  characterized  as  a 
flattened  density  distribution  that  dominates  locally  for  stars  with 
[Fe/H]  >  —  2,  whereas  the  outer-halo  population  is  nearly  spherical,  and 
dominates  at  larger  distances  and  locally  for  stars  with  [Fe/H]  <  —  2. 
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In  Table  1  of  this  Article,  rows  12  to  20  (the  ‘Field  F,  G  turnoff  (non-kinematic)’) 
columns.  The  corrected  table  is  shown  below. 

Table  1  |  Studies  claiming  a  retrograde  outer  halo 

were  inadvertently  moved  up  one  row  in 

the  N  and  (V^) 

Sample  and  selection  criteria 

N 

Additional 

restrictions 

(kms  1) 

Method 

Source 

Globular  clusters  (non-kinematic) 

19 

'Young  halo' 

-64  ±74 

F&W 

Ref.  2 

Globular  clusters  (non-kinematic) 

20 

'Young  halo' 

-42  ±80 

F&W 

Ref.  10 

RR  Lyrae  stars  (non-kinematic) 

26 

Z  <  8  kpc 

-95  ±  29 

FSM 

Ref.  9 

Field  subdwarfs  (kinematic) 

30 

Zmax  '>  5  kpC 

-45  ±22 

FSM 

Ref.  7 

Bias  corrected 

±24  ±13 

Ref.  16 

Field  horizontal-branch  stars  (non-kinematic) 

90 

[Fe/H]  <  —1.6 

-93  ±  36 

F&W 

Ref.  8 

|  Z  >  4  kpc 

Field  subdwarfs  (kinematic) 

101 

V<  -100  km  s_1 

-32  ±  10 

FSM 

Ref.  13 

[Fe/H]  <-1.8 

Field  F,  G,  K  dwarfs  (non-kinematic) 

250 

|Z|  >5  kpc 

-55  ±  16 

FSM 

Ref.  6 

Field  F,  G  turnoff  (non-kinematic) 

Zmax  >  5  kpc 

FSM 

This  work 

2,228 

[Fe/H]  <-1.0 

-11  ±2 

200 

[Fe/H]  <-2.2 

-41  ±  11 

Zmax  >  10  kpc 

771 

[Fe/H]  <-1.0 

-38  ±5 

94 

[Fe/H]  <  -2.2 

-71  ±  17 

^max  •'>  15  kpc 

371 

[Fe/H]  <-1.0 

-56  ±8 

54 

[Fe/H]  <-2.2 

-71  ±25 

216 
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